Resonant enhancement of ultracold photoassociation 
rate by electric field induced anisotropic interaction 

Debashree Chakraborty^ , Jisha Hazra^ and Bimalendu Deb^'^ 

^Department of Materials Science, Indian Association for the Cultivation of Science, 
Jadavpur, Kolkata 700032, INDIA 

^Raman Center for Atomic, Molecular and Optical Sciences, Indian Association for 
the Cultivation of Science, Jadavpur, Kolkata 700032, INDIA 

Abstract. We study the effects of a static electric field on the photoassociation of 
a heteronuclear atom-pair into a polar molecule. The interaction of permanent dipole 
moment with a static electric field largely affects the ground state continuum wave 
function of the atom-pair at short separations where photoassociation transitions occur 
according to Franck-Condon principle. Electric field induced anisotropic interaction 
between two heteronuclear ground state atoms leads to scattering resonances at 
some specific electric fields. Near such resonances the amplitude of scattering wave 
function at short separation increases by several orders of magnitude. As a result, 
photoaasociation rate is enhanced by several orders of magnitude near the resonances. 
We discuss in detail electric field modified atom-atom scattering properties and 
resonances. We calculate photoassociation rate that shows giant enhancement due 
to electric field tunable anisotropic resonances. We present selected results among 
which particularly important are the excitations of higher rotational levels in ultracold 
photoassociation due to electric field tunable resonances. 
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1. Introduction 

Research interest in the field of cold molecules has witnessed a tremendous growth in 
recent times. The possibility to investigate molecular behaviour at low temperatures 
motivates physicists and chemists from diverse backgrounds to study cold and dense 
molecular gases. Certainly, this development is inspired by the great success in the 
closely related field of cold atoms. Molecules can have properties which are not available 
with atoms, for instance a hetero-nuclear molecule can possess permanent electric dipole 
moment. While the rich internal structure of molecules possesses new challenges for 
cooling and trapping them, the potential applications of cold molecules are remarkable. 
This is especially true for ultracold polar molecules which allow for a wealth of interesting 
studies. The interaction of permanent dipoles with an external electric field provides a 
tool to manipulate physical processes taking place in the ultracold regime. Several 
quantum computation devices have been proposed based on the interaction of the 
heteronuclear dimers with external electric fields [U [2]. External electric fields can 
also alter the internal rovibrational structure and dynamics of heteronuclear molecules 
[21 m |5]. The orientation, angular motion, hybridization, and changes in the transition 
rates and lifetimes of the rovibrationally excited state of a LiCs molecule in a strong 
static electric field have been theoretically studied |6l |7]. Electric field is used to 
enhance the interaction between Li and Cs atoms in an ultracold collision [HI [9]. 
This enhancement is due to the interaction of the instantaneous dipole moment of 
the heteronuclear collision complex with the static electric field. The interaction of a 
heteronuclear pair of atoms with an electric field can couple states of different angular 
momenta in the electronic ground continuum of the pair leading to resonances in an 
ultracold collision. 

Our purpose here is to utilize the static electric field induced scattering resonances 
to infiuence free-bound photoassociation (PA) [lOj process of molecule formation. PA of 
ultracold atoms via the interaction with an electromagnetic field has become a standard 
technique to produce cold and ultracold molecules. The experimental techniques 
previously used for the formation of homonuclear molecules have been extended to 
heteronuclear alkali-metal dimers. Several experimental groups have reported the 
photoassociative formation of ultracold alkali-metal dimers , such as NaCs [Tl], KRb 
[121 Ha [H], RbCs [la US] and LiCS [T7] in their electronic ground states. The effect 
of a static electric field on the formation of heteronuclear molecules in their electronic 
ground state via one-photon stimulated emmision from ground continuum has been 
theoretically studied [T8] . 

Here we investigate the effect of a static electric field on the photoassociation process 
of a heteronuclear atom pair to produce molecules in an excited electronic state. Since an 
electric field can couple different angular momentum states (partial waves), we need to 
investigate anisotropic scattering at low energy. Taking LiCs molecule as a prototype, 
we first present a detailed investigation of low-energy anisotropic atomic collision in 
the presence of a static electric field. We have used Numerov-Cooley algorithm-based 
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Figure 1. Adiabatic Born-Oppenheimer potentials of LiCs without hyperfine 
interaction for excited (a) and ground (b) states. The ground state asymptotically 
corresponds to both Li and Cs atoms in electronic S state while excited state 
corresponds to Li in S state and Cs in P state. 



multichannel scattering techniques and found several anisotropic resonant structures in 
the scattering cross-sections of ^Li + ^^^Cs collision. Although the effects of a magnetic 
Feshbach resonance [19] on PA has been recently studied both experimentally [201 121] 
and theoretically [221 l23l [211 |25l |26] , to the best of our knowledge, the effects of electric 
field induced anisotropic resonances on PA into excited molecular levels have not been 
yet studied. One notable feature of the influence of these anisotropic resonances on 
PA is the occurrence of higher rotational excitations in excited molecule. This follows 
from the large modification of two-atom continuum wave functions for higher partial 
waves. Owing to the strong spatial dependence of permanent dipole moment at short 
separations, the interaction of a static electric field with the permanent dipole moment 
of a heteronuclear atom-pair leads to large enhancement of continuum wave functions 
at short separations, near electric fields at which anisotropic resonances occur. 

This paper is organised as follows: First, we briefly discuss how an external 
static electric field modifies the effective interaction potential between two heteronuclear 
collision pair. In section 2, we present mathemetical formulation of the problem with 
an emphasis on anisotropic scattering. The effects of an external static electric field on 
the ground state scattering of "^Li + ^^^Cs are discussed in section 3. In section 4, we 
discuss the effects of static electric field-induced anisotropic interaction on PA rate and 
present our main results. Finally we summarize and come to the conclusion in section 
5. 
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Figure 2. The potentials V(R) of LiCs in different diabatic hyperfine channels as a 
function of internuclear distance R are plotted. Panel (a) shows the short range and 
panel (b) shows the long range part of the potentials, /i and /2 are the hyperfine 
quantum numbers of ^Li and ^'^'^Cs, respectively (See Tabel-I in the text). 



2. Formulation 

The dynamics of Li-Cs collision in the presence of an electric field is effectively described 
by the radial Hamiltonian 

where /i = 17111712/ (mi + 7712) is the reduced mass of two atoms ^Li and ^'^^Cs with masses 
rui and m2 respectively, R is the interatomic distance, is the rotational angular 
momentum of the collision complex and the angles 6 and specify the orientation of 
the interatomic axis in the space fixed coordinate frame. The electronic interaction 
potential can be represented as 

VciR) = EE I SMs)VsiR){SMs I (2) 

5 Ms 

where S is the total electronic spin of the two atoms and Ms is the projection of S on 
the Z-axis. Vs{R) represents the adiabatic interaction potential of the molecule in the 
spin state S. The ground state Vg{R) and the excited state VexiR) potentials as shown in 
figure 1, are taken from [27] and [28], respectively. This interaction is therefore diagonal 
in the adiabatic basis | IMj; SMs), 

{S'M's; I'M', I V,{R) \ SMg; IMj) = 5ij,5MjM',^s,s'hisM'sMR) (3) 

where S = si + S2 and I = ii + i2, si and S2 are the electronic spins while ii and 12 are 
the nuclear spins of the Li and Cs atoms, respectively. The total hyperfine Hamiltonian 
Vhf for the collision complex can be written as a sum of the two atomic hyperfine 
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Figure 3. Total elastic cross-section a in a.u. vs electric field £ in kV/cm is plotted for 
different total number of partial waves to show the convergence of the first resonance 
point at f = 1298 kV/cm. 
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% = E ^^l-^- = E ^ifl - ^1 - (4) 

j=i n zn 

where fj = Sj + ij is the total spin and a^flj is the hyperfine constant of atom j which 
is 402.00 MHz for ^Li and 2298.25 MHz for ^^^Cs |29]. This hyperfine interaction is 
diagonal in the atomic or diabatic or long range basis | i, f, fi, f^) in which K(-R) can 
be expressed as, 

= E ys{SMs; IMi- imt I iM2) fmfime) 

S' Ms,, I'M J, 

X ((/(/2)/ V; ^'^i' I S'Msr, I' Mir, £'m,')- (5) 
The transformation from coupled hyperfine representation to short range representation 
is given by 



6u'Sm,mASMs;IMj I /m^)^(2/i + l)(2/2 + l)(25+l)(2/+l) 

51 ii fi 

52 12 f2 ? (6) 

S I f _ 

where (SMs; IMj \ frrif) is Clebsch Gordon coefficient and the quantity in curly braket 
is 9j-symbol. The hyperfine Hamiltonian can only couples channels with the same total 
angular momentum projection ruf = nifi + mf2 = Ms + Mj where rrifi, mf2, Ms 
and Mj are the projections of /i, /2, S and I respectively. For ''Li, ii = 3/2 and for 
""^^^Cs, ^2 = 7/2 (both with si = S2 = 1/2). Since rrif takes values from -6 to +6, 
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Table 1. Three different hypcrfine channels of ^Li and ^^'^Cs. 



Channel index 


/i('Li) 




/ 


a 


2 


4 


5 


b 


1 


4 


5 


c 


2 


3 


5 



the total degeneracy of the atom pair is 128. However, atoms mainly collide on the 
Li(2^S'i/2, /i = 2) + Cs(6^S'i/2, /2 = 3) channel, and thus only 35 degenerate entrance 
channels have to be considered. For / = 5, three coUisional channels are possible as given 
in Table 1. The potentials of these three channels are shown in figure 2 . We consider 
a single asymptotic hyperfine channel (c) with /i = 2, /2 = 3 since the difference in 
energy of this channel from the other two channels is large enough compared to collision 
energy at ultracold temperatures. Thus we can approximate our calculation as a single 
channel with fi = 2, /2 = 3. The operator Vs^R) describes the interaction of the atoms 
with an external electric field. It can be written in the form 

VsiR) = -ScosOY^Yl I SMs)ds(R)(SMs | (7) 

S Ms 

where d5(R) denotes the dipole moment functions of LiCs in the different spin states 
and S the electric field magnitude. Li and Krems [9] have given an analytical expression 
for this dipole moment function approximating the numerical data computed by Aymar 
and Dulieu [30] . This analytical expression is given by 

dsiR) = D exp [-a{R - Ref] , (8) 

with the parameters Re = 7.7 oq, a = 0.1 Oq^, and D = 6 Debye for the singlet state 
and D = 0.5 Debye for the triplet state (where Oq = Bohr radius ). The matrix element 
of V£{R) are evaluated using the expressions 

{em, I COS ^ I i'm,,) = 5m,m,{-ir^ ^ L n n (9) 

\—mi niii / \0 0/ 

X [(2£+l)(2£' + l)]'/' 

and 

{SMs I (eE I S"M';,)ds»{S"M'^ I ] I S'M's) = dsdss'hisM^, (10) 

\S" M'^ J 

where {R \ i,mi) = YimiiR)- Due to the symmetry of the z-component of the angular 
momentum, the matrix element {imi \ cos 6 \ i'mgi) exists only if £ — £' = ±1 and 
me = mil. Therefore electric field can couple even and odd parity channels to each 
other but not to themselves. 

The scattering wave function can be expressed as 

^(^) = ^E^i,n^,iR)Y;„.,ik)YernAR) (H) 
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Figure 4. Energy-normalised scattering wave function (/)o,o (a.u.) is plotted as a 
function of interatomic distance R (a.u.) at the first resonant electric field £ = 1298 
kV/cm (a), £ = 1200 kV/cm (b) and £ = (c). 



which has the asymptotic form 

m{R) ~ Rexp{ik.R) + fj^p exp(iA;i?), (12) 
where k and k' are the incident and the scattered momentum, respectively. The on- 

— * ^ 

shell elastic scattering is then described by /^r^,, with the scattered momentum k' — kR. 
Expanding the scattering amplitude /^p into the complete basis 1^^^, we have 

k,k'-^T.T,^cik)Y,^e{R) (13) 
where T^^^ is a T-matrix element which can also be expanded in the following form 

^-Ti^x^) = E (14) 



and 



^ 01 



eMik.R) = 47r i'M{kR)Y;^fi)Y,,^^{R). (15) 



Substituting equations (13) and (15) in equation (12) we get 

471 



R-^iR) ^^Y: AYLfi) sm{kR - in/2) + T.^^^) 



X e-K^{ikR-ilT:/2)]Ye,mXR) (16) 
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Figure 5. Plot of stimulated linewidth F vs. electric field £ at laser intensity IW/cm^. 



where the asymptotic form jt{kr) ~ sm.{kR — i7r/2)/{kR) has been used. We get the 
multichannel form of the scattering equation which is given by 



hejPe^,{R) = Yl i''~H^rm \ Vs{R) \ f m,.)# 



where 



he — 



+ 



+ 1) 



+ V -E. 



(17) 



(18) 



E = h'^k^/2ij, is the collision energy and V is the central potential for the chosen hyperfine 
channel. The asymptotic boundary condition on ipemi can be set as 

^em, ~ sm{kR - £7r/2) + Te„,, (k) cxp{tkr - Ut:/2). (19) 

Under the approximation of single asymptotic hyperfine channel as discussed earlier, 
we can write the coupled equation (17) in a matrix form in relative angular momentum 
(£) basis as given by 



(20) 



where I is the identity matrix, = Diag[l(l+1), 2(2+1), , £{£+!) ] is a diagonal 

matrix. The wave function $ is a matrix whose elements are given by (f)Y-y where 
•J = i,me and 7' = f, me' are the incident and the scattered angular states respectively. 
Therefore the overall wave function for an incident partial wave 7(£, m^) becomes 

^^Jy=e,me{E, R)^ E R) I t')- (21) 

By imposing the boundary conditions on the partial waves as given by equation (19), 



we can get the T-matrix elements, f 



. The total elastic cross-section is given by 



(22) 
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Figure 6. Square of franck-condon factor | rj p (a.u.) vs. electric field £ (kV/cm) 
for transitions between the excited bound state v = 26, J = 1 and different scattering 
states (pi'oio with incident partial wave £ = and scattered partial wave i' = (solid 
line), i' = 1 (dashed line) and £' = 2 (dotted line). 

3. Ground state scattering: Electric field induced resonances 

To obtain asymptotic scattering solutions, Numerov technique was adopted to 
numerically propagate equation (20). In the asymptotic region, the free wave function 
goes hke 

(f)i^ = ci sin(A;i? - In/ 2) + C2C0s{kR - Itt/2). (23) 

We construct 5- matrix with its elements given by (ci —ic2)~^{ci+ic2) and the T matrix 
{S — l)/(2i). In the absence of electric fields, different partial wave states of the Li-Cs 
collision complex are uncoupled and s-wave scattering almost entirely determines the 
cross-sections at ultralow kinetic energies. The interaction of pemanent dipole moment 
of the colliding pair of heteronuclear atoms with electric field as given by equation (7), 
induces coupling between different angular momentum states and may thus affect the 
ground state scattering wave functions. Ultracold s-wave scattering is isotropic: the 
probability to find the atoms after s-wave collisions does not depend on the scattering 
angle. The interaction with electric fields, however couples the spherically symmetric 
s-waves to anisotropic p-wave which in turn is coupled to d-wave, d-wave to g-wave and 
so on. Figure 3 shows anisotropic resonance structures in the scattering cross-section 
of ^Li -I- ^^^Cs collision at 50//K energy. Typically convergent results are obtained for 
a minimum angular momentum of ^= 9, as shown in figure 3. Larger the electric field, 
larger is the number of partial waves required for convergence. We get the first resonance 
peak near 1298 kV/cm and the second one near 1650 kV/cm. To show the effect of these 
resonances on PA, we have carried out PA calculations near the first resonance peak. 
The permanent dipole moment function of the collision complex as given by 
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Figure 7. The photoassociation rate KpA (in unit of 10 ^^cm'^s ^) is plotted as a 
function of atom-field detuning d in GHz £ = (a) and at resonant electric field £ = 
1298 kV/cm. 



equation (8) is typically peaked around the equilibrium distance of the diatomic molecule 
in the ground state and quickly decreases as the atoms separate. Due to this strong 
spatial dependence of permanent dipole moment at short separations, the interaction of 
a static electric field with the permanent dipole moment of the collision pair leads 
to large modification of continuum wave functions at short separations. Figure 4 
demonstrates that the modification is particularly significant near the resonant electric 
fields. Since PA transitions take place at short separations near the outer turning point 
of B^U, modification of the continuum wave functions at short separations influences 
PA transition probability as described in the following section 4. 



4. Photoassociation: Electric field effects 



Recently, an experiment has demonstrated the formation of ultracold bosonic ^Li^^^Cs 
molecule in their rovibrational ground state by a two-step PA procedure [31]. For 
illustration of electric field effects on PA, we consider PA transition from Li(2^S'i/2, /i = 
2) + Cs(6^S'i/2, /2 = 3) continuum to the f = 26 , J = 1 level of B^U, near the electric 
fields at which anisotropic resonances occur. Photoassociation rate coefficient [32] is 
given by 

Kpa{T,ul) = ^^£(2^+1) I SpA{E,e,coL) (24) 

where Vrei = hk/ fi is the relative velocity of the two atoms, Spa{E, i, up) is the S-matrix 
element for the process of loss of atoms due to PA and (....) implies an averaging over 
thermal velocity distribution. Assuming Maxwell-Boltzman distribution at temperature 
T we have, 

Kpa{T,ul) = -—Y,{2i+l) / I SpA{E,i,UL) r e(-^/^^^)c^E (25) 
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Figure 8. The photoassociation rate Kpa (in unit of 10~^''cni^s~^) is plotted as a 
function of atom-field detuning 6 in GHz at the resonant electric field £ = 1298 kV/cm. 
The rate increases for J = 3, 4, 5 and 6 rotational levels and then start decreasing for 
J> 7. 



where 

I SpA{E,e,ujL) |'= p^-p — . (26) 

[(5 + ^)^ + ((r + 7s)/2)2] 

Here Qt = {2'kijlKbT /h^Y^"^ is the translational partition function and 7^ is the natural 
linewidth of the photoassociated level. E = is the relative kinetic energy of 

the colliding pair of atoms with reduced mass 11, Ey^j is the bound state energy and 
5 — ujl — oja '^s the frequency offset between the laser frequency ujl and atomic resonance 
frequency uja- The stimulated linewidth F is given by 

r = — I I Dt{R) I il^^=i,mi{E, R)) |2 (27) 

where I is the intensity of the PA laser, eo and c are the vaccum permitivity and speed 
of hght, respectively. Dt(R) is the transition dipole moment. At zero electric field and 
for 50;LiK collisional energy we find F = 0.315 kHz at laser intensity IW/cm^. Figure 
5 demonstrates that the stimulated linewidth as a function of electric field S has a 
resonant structure with giant enhancement near the resonant electric field even at low 
laser intensity. This is because of large modifications of continuum wave functions at 
short separations. To know which components (f)^^^ of the continuum wave functions 
have significant contribution to such resonant enhancement of F, we plot square of 
Pranck-Condon (FC) overlap integral | 77 p for a few wave component (j)e'o,m in figure 6. 
77 is defined by 

77 = {(f>v,j I ^^^^=e,m^{E, R)) = | ^ R) \ 7'). (28) 

'y'=l',m[ 

In fact, we find that a large number of partial waves significantly contribute to the 
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Figure 9. (a) The unit normalised bound state (py^j in a.u. In panel (b) plotted 
are the energy normalised scattering wave functions 0^',^=o as a function of R (a.u.) 
for e = Q (dashed-dotted line), I' = I (dotted Une), 1' = 2 (sohd line) and ^' = 3 
(dashed line). Panel (c) shows (f)e'^e=o for £' = A (dotted line), £' = 5 (solid line), £' = 
6 (dashed-dotted line) and i' = 7 (dashed line) at the resonant electric field £ = 1298 
kV/cm. 



resonant features. This is perhaps due to the formation of a quasibound complex inside 
the centrifugal barrier of a large number of partial waves which are strongly coupled 
among them. Figure 7 shows that KpA near the resonant electric field increases by four 
orders of magnitude compared to that at zero electric field. This enhancement is due 
to the electric field induced coupling between different partial waves. The coupling is 
so strong that higher rotational levels of the excited electronic state can be populated 
as shown in figure 8. These higher rotational excitations are not possible at ultracold 
temperatures in the absence of electric field. The PA rate increases for rotational levels 
J = 3, 4, 5 and 6 and then starts decreasing for J > 7 as shown in figure 8. This can 
be understood from figure 9 which shows that the amplitude of (j)e'=4/=o, 0^'=5/=o and 
0£/=g £=0;Wave functions near the outer turning point of the excited molecular level are 
larger than 0^/=i,£=o, (j)t=2,e=o and 0£/=3£=o waves. The selection rule for free-bound 
transition J — \ L + S\\ <£' < |j+ |^ + ^|| allows the scattered partial waves £' — 
0, 1 and 2 to be accessible for PA with J=l. Similarly, higher rotational levels J > 1 
become populated due to excitations of higher partial waves in the ground continuum 
through electric field induced anisotropic scattering. 

The enhancement of PA rate near resonant electric fields critically depends on the 
outer turning point of the molecular bound state that is accessible to the PA transition. 
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As mentioned earlier, the anisotropic resonances bring about large modification of 
the amplitude of continuum wave functions at relatively short separations. Now, if a 
prominent antinode of the modified continuum wave functions lies at a separation near 
the outer turning point of the bound state, then enhancement of PA rate is expected due 
to large FC overlap. We have chosen a particular bound state for which enhancement 
can be achieved. Although, in this paper we have studied enhancement of the rate of 
formation of excited molecular states only, the similar method of electric field induced 
resonances can readily be extended to enhance the formation rate of ground molecular 
states via one photon stimulated emission [18]. Enhancement can also be achieved 
using quantum interference [221 ESI |33] via laser-coupling another bound state with 
the continuum. In case of PA in the presence of a magnetic Feshbach resonance, a 
quasi-bound state embedded in the ground continuum gives rise to Fano type quantum 
interference p3] leading to enhancement in PA pO]. In case of heteronuclear atoms, 
quantum interference can be used to enhance the production of ground molecules due to 
the existence of permanent dipole moment. Recently, an all optical method of quantum 
interference scheme has been proposed for efficient production of ground state polar 
molecules [33]. The unique feature of the electric field induced enhancement would be 
the controllability of the excitations of higher rotational states. 

5. Conclusion: 

In conclusion, we have analysed the effects of a static electric field on the PA of a 
heteronuclear atom-pair. Our results show that it is possible to enhance PA rate 
by several orders of magnitude by tuning electric field near anisotropic resonances. 
Due to anisotropic nature of the interaction, a large number of partial waves in the 
ground continuum become strongly coupled. As a consequence, higher rotational levels 
can be populated in an excited dimer formed by PA at ultralow tempertures. This 
leads to the possibility of selectively populating higher rotational levels in ground 
state polar molecule by Raman-type two-color coherent PA. Furthermore it may be 
interesting to investigate into the effects of electric field on stimulated Raman adiabatic 
passage (STIRAP) from free atoms to ground state molecules. Electric field induced 
resonances may be coupled with multiple transition pathways to devise novel opto- 
electrical quantum interference schemes for efficient production of selective rovibrational 
molecular states. 
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